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Abstract 

We discuss the effects of fast nuclear excitations on heavy-ion fusion reactions at 
energies near and below the Coulomb barrier. Using the fusion of two ^°Ca nuclei as an 
example and the inversion method, we show that the mass renormalization induced by 
fast nuclear excitations leads to a large surface diffuseness in the effective potential for 
heavy-ion fusion reactions. 
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Heavy-ion fusion reactions at energies below the Coulomb barrier are typical examples 
of the macroscopic quantum tunneling, which has been a very popular subject in the past 
decade in many helds of physics and chemistry [jl|, One of the major interests in the 

macroscopic quantum tunneling is to assess the effects of environments on the tunneling 
rate of a macroscopic variable. When the bombarding energy is below the Coulomb 
barrier, heavy-ion fusion reactions take place by a quantum tunneling. It is now well 
established that the fusion cross section at sub-barrier energies is enhanced by several 
orders of magnitude compared with the prediction of a one-dimensional potential model 
due to the coupling of the relative motion to nuclear internal degrees of freedom, which 
play the role of environments 0. The present authors have shown that the coupling of the 
tunneling degree of freedom to fast environmental degrees of freedom causes both a static 
potential and a mass renormalizationsj^. On the other hand, the hrst and the second 
derivatives of the experimental fusion cross section with respect to the bombarding energy 
suggests that the fusion reaction between two ^°Ca nuclei can be described by a potential 
model using the surface diffuseness parameter a=1.27 fm^, 0, which is much larger than 
the usually accepted value a ~ 0.6 fm. The spin distribution of the compound nucleus 
made by the fusion of with also suggests a large effective surface diffuseness 

parameter a=1.27 fm[Q. Such a large surface diffuseness has already been suggested long 
time ago through the inversion method of obtaining the effective fusion potential directly 
from the experimental data. In that procedure, a much larger surface diffuseness a=1.59 
fm was suggested for ^°Ca -|- ^°Ca scattering 0. 

The aim of this paper is to discuss the connection between this large effective surface 
diffuseness and the potential and the mass renormalizations induced by the coupling of 
the relative motion between heavy ions with fast nuclear intrinsic degrees of freedom. We 
consider the fusion of two ^°Ca nuclei as an example. We hrst determine the bare nuclear 
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potential in the entrance channel by the double folding procedure. We then calculate the 
excitation function of the fusion cross section by taking the effects of a few vibrational 
modes of excitation of ^°Ca into account through the adiabatic potential as well as the 
mass renormalizations. We then use the inversion procedure to obtain the effective fusion 
potential from thus obtained excitation function of the fusion cross section and discuss 
its surface properties. 

In the double folding procedure, one often expresses the densities of nuclei by a step 
function of range Rd convoluted with a Yukawa smearing function of range l/^d and 
assume a Yukawa interaction for the nucleon-nucleon interaction with the range parameter 
K/. The resulting potential consists of two parts, which are proportional to and 

q-kiR surface region[^. If the Michigan three range Yukawa (M3Y) interaction |[TT|, 


is assumed as is often done for the nucleon-nucleon interaction, kj is larger than 


the standard value of Kd- The diffuseness parameter of the double folding potential is 
then approximately given by the inverse of Kd- More precisely, the surface diffuseness 
parameter slightly deviates from this simple consideration. Actually, Akyiiz and Winther 
numerically performed the double folding procedure and obtained the value 1/1.36 fm 
for the surface diffuseness parameter in the potential between two ^°Ca nuclei when they 
assumed the charge distribution of ^°Ca to be given by ^^=1-55 fm“^ and used the M3Y 
force as the nucleon-nucleon interaction]^. Note that 1/1.36=0.75 fm is slightly larger 
than the standards of what is usually used in optical-model calculations, but not so much 
as those suggested from the analyses of the experimental sub-barrier fusion cross section. 
Henceforth we take the Woods-Saxon shape for the bare nuclear potential between two 
"^^Ca nuclei 


VNiR) = 


( 1 ) 


_j_ ^{R ^Ro)/^ 

and use the value 0.75 fm for the surface diffuseness parameter a as suggested in the above 
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arguments. Following ref.[|^, the strength and the range of the bare nuclear potential 


are taken to be 


Vo = IQ'K'yaR 

Ro = 1 . 2 - 40 ^/ 3 - 0 . 09 , 


( 2 ) 

(3) 


where 7 and R are 0.95 and Ro/2, respectively. The bare nuclear potential thus obtained 
gives the Colulomb barrier at lO.Ofm with the height and the curviture 53.1 MeV and 
3.43 MeV, respectively. 

We now consider the effects of surface vibrations of '^^Ca during the fusion process. 
We ignore the effects of nucleon transfer reactions, althogh they play an important role 
in quantitatively explaining the experimental fusion cross section |l^, Following 

ref.[|^ we consider the excited states of ^°Ca at 3.74, 3.91, and 4.49 MeV with spin- 
parity 3“, 2~^ and 5“ to be the one phonon state of three different vibrational modes of 
excitation, though the 2^^ state will rather be the 2^^ member of the rotational band based 
on the O’*" state at 3.35 MeV[^. Our interest is to study the effects of the coupling of 
the relative motion to these states on the surface diffuseness parameter in the effective 
potential for fusion reactions. 

We assume the linear oscillator coupling and use the no-Colioris approximation in 


order to reduce the dimension of the coupled-channels equations[0, |I^]. The effective 
Hamiltonian for this system then reads 


H = 






+ 


2 d \ J{J + l)h^ 


2M \dR^ RdR 2MR^ 


+ Vc{R) + Vn{R) 


E Uwa (+ ^ ) + 

A=3,2,5 ' ^ 


/2A + 1 


dvr 


(do'/jv(-R) + /c(fl)) (7 + oa) (4) 


where M, R, J, VciVj^, and A are the reduced mass, the relative distance, the total angular 
momentum, the Coulomb and the nuclear interactions between the two ^°Ca nuclei, and 
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the multipolalities of the vibrational excitations, respectively. fN{R) and fc{R) are the 
nnclear and the Coulomb coupling form factors, respectively. We determine them fol¬ 
lowing the collective model|^. hu}\, a^iax), and are the excitation energy, the 
creation (annihilation) operators of the oscillator quanta, the nuclear and the Coulomb 
coupling strengthes of the 2^-pole vibrational excitations, respectively. The factor 
in front of the coupling strength arises from the no-Coriolis approximation. We take the 
coupling strengthes from ref.[|^ (see Table I). In Table I, Rc and Rn are the charge and 
the matter radii of “^^Ca, respectively. The Coulomb coupling strength was determined 
from the transition strength B{EX) from the one phonon state of each vibrational mode 
of excitation to the ground state. The nuclear coupling strength was modihed to repro¬ 
duce the data of inelastic scattering of from ^°Ca. All the coupling strengthes are 
multipled by a/2 in order to take the mutual excitations of the target and the projectile 
into account. 

The excitation energy huj of all the vibrational excitations which we are interested 
in now is larger than the curvature of the bare potential barrier. Therefore, the present 
problem corresponds to the case of an adiabatic tunnelig, i.e. a slow tunneling. The 
adiabaticity of the tunneling process depends on two parameters [pO|] . One is the ratio 
of the energy scale Oi = hfl/hu, hfl being the curvature of the bare potential barrier. 
The other is the ratio of the strength of the coupling Hamiltonian to the excitation 
energy 02 = aocRa/hoj, Ra and aoc being the thickness of the tunneling region and the 
derivative of the coupling Hamiltonian in the case when one assumes a bi-linear coupling, 
respectively. The adiabatic assumption breaks down if these parameters are larger than 
one. The typical values of these parameters for ^°Ca -|- ^°Ca scattering are listed in Table 
H. In order to obtain 02 , we estimated the derivative of the coupling Hamiltonian at the 
position of the bare potential barrier and determined Ra at energy 1 MeV below the top 
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of the bare potential barrier based on the parabolic approximation. Table II justihes the 
use of the adiabatic formula to calculate the barrier penetrability for the fusion reaction 
in the present system. We use the revised adiabatic formula of ref.[^]. It modihes the 
well known adiabatic formula, which takes only the static potential shift into account, by 
adding the mass renormalization. Our formula for the barrier penetrability then reads 


with 


S{E) = {VUR) - E) 


( 6 ) 


for the s-wave scattering. In eq.(6) Ri and R 2 are the inner and the outer classical turning 
points, respectively. The renormalized potential and mass in this equation are given by 

§. 


Vm(R) = Vc{H) + Vk{R)- Y. 


A=3,2,5 


2A + 1 
47r 




2 1 




2A + 1 

Mad{R) = M + 2 — -l« 


A=3,2,5 


dvr 


iN) dfN{R) , iC) dfc{R) \' 
AO 7 D ' *^A0 7 D I 


(7) 

( 8 ) 


dR ' dR J hul 
Fig.l shows the tunneling probability in the s-wave scattering as a function of the 
bombarding energy. The dotted line is the penetrability when there is no coupling, i.e. 
the penetrability for the bare potential barrier with the bare mass M. The dashed line is 
obtained with the standard adiabatic formula, where only the adiabatic potential renor¬ 
malization is taken into account. The result of the revised adiabatic formula which takes 
the mass renormalization into account is denoted by the solid line. Notice that the overen¬ 
hancement of the effects of the channel coupling when one considers only the adiabatic 
potential shift (the dashed line) is corrected by the effects of the mass renormalization. 
Fig.2 shows the renormalized mass Mad{R) in this system (eq.(8)) in the ratio to the bare 
mass M. It is a function of the separation distance between the colliding nuclei reflecting 
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the radial dependence of the coupling form factor. It has two peaks because of the strong 
cancellation between the Coulomb and the nuclear couplings. The small spike at around 
i?=8fm originates from the discontinuity of the derivative of the Coulomb coupling form 
factor. This defect plays, however, no essential role. 

We now derive the effective fusion potential corresponding to the solid line in Fig.l. 
We thus demonstrate how the effects of the channel coupling including the mass renor¬ 
malization can be mocked up by renormalizing the surface diffuseness. To this end, we 
use the inversion procedure which gives the thickness of the tunneling region, i.e. the 
thickness of the potential barrier, at each energy directly from the excitation function of 
the fusion cross section||] 


m = -- 




rB 


dS{E') 

-.dE' 


(9) 


2M Je ^JE' - E 

where S{E) is the action integral given by eq.(6) and B is the barrier height detemined 
from the condition S{B) = 0. This is the height of the potential barrier including the 
adiabatic potential shift. Note that the mass M in the factor in front of the integral in 
eq.(9) does not include the mass renormalization, but the bare mass. This is because our 
aim is to get the effective potential for the relative motion with bare mass. The results 
of the inversion procedure are shown in £g.3 by the solid line. For the integration in 


eq.(9), we changed the variable from E' to y/E' — E to eliminate the singularity at the 
lower limit of the integral]^. The inversion procedure gives only the barrier thickness 
t{E). Consequently, the turning points themselves as functions of the energy remain 
indeterminate. We assumed that the outer turning points are given by those of the 
adiabatic potential barrier and determine the inner turning points of the effective potential 
barrier by using the barrier thickness given by the inversion procedure. For comparison, 
£g.3 also shows the bare potential barrier (the dotted line) and the adiabatic potential 
barrier (the dashed line). 


7 








Since we are interested in the surface diffuseness parameter of the nuclear potential, 
we subtract the Coulomb potential Vc{R) from both the effective potential obtained by 
the inversion procedure and the adiabatic potential. We call these potentials the effective 
nuclear potentials, because they contain not only the effects of the nuclear coupling but 
also those of the Coulomb coupling. Fig.4 shows the effective nuclear potentials obtained 
in this way. The dashed line denotes the effective nuclear potential in the adiabatic limit, 
i.e., the potential obtained by subtracting the Coulomb potential from the adiabatic 
potential. The effective nuclear potential, which takes account of the mass renormaliztion 
as well, is denoted by the solid line. Also shown are the potentials of an exponential shape 
with two different surface diffuseness parameters. Their strength was chosen, such that 
the values of the resultant potentials agree with that of the adiabatic potential at the 
peak position. The dot-dashed and the dotted lines correspond to the potential when the 
surface diffuseness parameter was chosen to be 0.7 and 0.8 fm, respectively. These choices 
simulate the effective potential in the adiabatic limit (the dashed line), and that which 
includes the effects of the mass renormalization(sohd line). Remember that the surface 
diffuseness parameter for the bare potential is about 0.75 fm. Fig.4 shows that the surface 
diffuseness parameter is smaller than that of the bare nuclear potential if only the effects 
of the adiabatic potential shift is taken into account, but gets larger than that of the bare 
nuclear potential when the effects of the mass renormalization is added. We thus conclude 
that one of the origins of the large surface diffuseness of the effective nuclear potential 
in heavy-ion fusion reactions is the effects of the coupling to nuclear inelastic excitations, 
which lead to both the potential and the mass renormalizations. 

In summary, we discussed the effects of nuclear intrinsic excitations on heavy-ion fusion 
reactions. We recast the potential and the mass renormalizations due to fast inelastic 
excitations of scattering nuclei in terms of the renormalized surface diffuseness in the 



effective fusion potential. As a concrete example, we analyzed the fusion reactions between 
two ^‘^Ca nuclei. We used the inversion method to obtain the effective nuclear potential 
and found that the surface diffuseness parameter becomes larger than that obtained by 
the double folding method if both the adiabatic potential shift and the adiabatic mass 
renormalization are taken into account. The value of the surface diffuseness parameter we 
obtained is about 0.8 fm compared to the standard value 0.6 fm. Although this value is 
not sufficiently large to explain the large surface diffuseness empirically obtained from the 
experimental data of sub-barrier fusion cross sections, we conclude that the large surface 
diffuseness obtained from data analyses partly represents the effects of channel coupling. 
It would be interesting to study how the effects of transfer reactions further modify the 
surface properties of the effective fusion potential. The large surface diffuseness parameter 
has also been suggested in the fusion reactions of on ^®^Sm, which is a fast tunneling 
rather than a slow tunneling discussed in this paper [§. It remains an open question 
how we should take into account the effects of the channel coupling in such a system to 
explain the large surface diffuseness. Contrary to the case of slow tunneling discussed 
in this paper, one has to treat a distributed potential barriers, or an energy dependent 
effective potential barrier |2T[|. 
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Society for the Promotion of Science for Young Scientists. This work was supported by 
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Aid for Scientihc Research on Priority Areas,Contract No.05243101,from the Japanese 
Ministry of Education, Science and Culture. 
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Table I 



hcuA(MeV) 


i?7va^o^(fm) 

3 - 

3.737 

0.6576 

0.4455 

2+ 

3.905 

0.1952 

0.1768 

5 - 

4.492 

0.4837 

0.2475 


Table II 


A^ 

Oi 

0-2 

3 - 

0.9177 

0.6176 

2+ 

0.8783 

0.1997 

5 - 

0.7635 

0.3668 


Table Captions 

Table I: Low lying vibrational states of ^'^Ca. and are the excitation 

energy of the 2^-pole vibrational state, and its Coulomb and nuclear coupling strengthes, 
respectively. Rc and Rn are the charge and the matter radii, respectively. 

Table II: Adiabaticity parameters for the tunneling process. Oi is the ratio of the energy 
scales of the tunneling and the environmental degrees of freedom hQ/hu, while 02 the 
ratio of the strengthes of the channel coupling Hamiltonian to the excitation energy of 
the environment a^cRa/fioj. 
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Figure Captions 


Fig.l : Excitation function of the s-wave barrier penetrability for ^°Ca + ^“^Ca scattering. 
The dotted line is the penetrability in the absence of the channel coupling. The dashed 
line was obtained by the standard adiabatic formula which takes only the effects of the 
adiabatic potential shift into account. The solid line was calculated with the revised 
adiabatic formula by adding the effects of the mass renormalization. 

Fig.2 : The renormalized mass (eq. (8)) as a function of the separation distance between 
the colliding nuclei in the unit of the bare mass. 

Fig.3 : Fusion barrier obtained by the inversion procedure. The dotted and the dashed 
lines are the bare and the adiabatic potential barriers, respectively. The solid line includes 
the effects of the mass renormalization. 

Fig.4 : Effective nuclear potential for the fusion of two ^°Ca nuclei. The solid line was 
obtained by subtracting the Coulomb potential Vc{R) from the effective potential barrier 
in the presence of the mass renormalization. The dashed line is the effective nuclear 
potential in the adiabatic limit disregarding the mass renormalization. The dott-dashed 
and the dotted lines are the exponential potentials with the surface diffuseness parameters 
0.7 and 0.8 fm, respectively. Their depth was chosen to reproduce the barrier height of 
the adiabatic potential. The surface diffuseness parameters a=0.7 fm and 0.8 fm were 
hxed to simulate the nuclear potential without (the dashed line) and with (the solid line) 
the mass renormalization. 
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